Half-magnetization plateau stabilized by structural distortion in the 
antiferromagnetic Heisenberg model on a pyrochlore lattice 
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Magnetization plateaus, visible as anomalies in magnetic susceptibility at low temperatures, are 
one of the hallmarks of frustrated magnetism. We show how an extremely robust half-magnetization 
plateau can arise from coupling between spin and lattice degrees of freedom in a pyrochlore antifer- 
romagnet, and develop a detailed symmetry of analysis of the simplest possible scenario for such a 
plateau state. The application of this theory to the spinel oxides CdCr204 and HgC^CU, where a 
robust half magnetization plateau has been observed, is discussed. 

PACS numbers: 75.10.-b, 75.10.Hk 75.80.+q 



Spinels, with chemical formula AB2O4, are ubiqui- 
tous among magnetic oxides. Notable examples of such 
materials include Fe304, a system exhibiting frustrated 
charge order and ferrimagnetism Q , the ci-electron heavy 
fermion compound L1V2O4 0, and the frustrated S — 
3/2 antiferromagnet ZnCr 2 4 0. In all of these com- 
pounds, the B-site ion is magnetic, and much of the 
beautiful strangeness seen in the behaviour of these com- 
pounds can be traced back to the fact that the B-site 
ions form an acutely frustrated pyrochlore lattice, built 
entirely of corner-sharing tetrahedra (Fig. 

The geometric frustration of the pyrochlore lattice is 
so great that both the classical (S = 00) and quantum 
(S = 1/2) antiferromagnetic (AF) Heisenberg models 
are believed to remain magnetically disordered down to 
T = 0,0. In real spinel oxides, however, the ground 
state degeneracy associated with the frustrated lattice 
geometry is usually lifted by a distortion of the lattice. 
ZnCr204, for example, undergoes a transition from a 
paramagnet with cubic symmetry to a Neel ordered phase 
with tetragonal symmetry at T = 12K. 

Another, more progressive, means of reducing the 
ground state degeneracy of a frustrated AF is to apply 
a magnetic field. Fields h greatly in excess of the ex- 
change coupling J between spins will remove magnetic 
frustration altogether and cause the system to become 
ferromagnetic. At intermediate fields h ~ J , frustrated 
AF's frequently undergo a succession of phase transi- 
tions, with associated anomalies in their magnetic sus- 
ceptibility. Where one particular state remains stable for 
a finite range of fields, a plateau is seen in the magne- 
tization curve M(h). Magnetization plateaus have been 
predicted to occur in both triangular lattice and Kagome 
AF's 0,0 EH- In P ure s pi n models, such plateaus oc- 
cur as an "order from disorder effect" where quantum or 
thermal fluctuations select one of many possible classical 
ground states F° r this reason they are usually very 



fragile, and relatively difficult to observe in experiment. 

In this letter we consider the interplay between mag- 
netic field, spin, and lattice degrees of freedom in a 
Heisenberg antiferromagnet on the pyrochlore lattice. 
Our main result is that the coupling of applied magnetic 
field to lattice distortion provides an extremely efficient 
mechanism for stabilizing a robust half-magnetization 
plateau, with exactly three up spins and one down spin 
in each tetrahedral sub-unit of the lattice. In principle, 
many such states may arise; we develop a detailed sym- 
metry analysis of the simplest case, in which all of the 
tetrahedra which go to make up the lattice distort in the 
same manner. 

Our analysis is of direct relevance to spinel oxides such 
as ZnCr204, where the A-site ion is non-magnetic, and 
the octahedrally co-ordinated B-site ion has an exactly 
half-filled ty, g shell of d-clectrons, giving rise to a spin 
S = 3/2 moment. And, indeed, just such a plateau has 
been observed in recent high field magnetization mea- 
surements on the closely related Cr spinels CdCr 2 04 and 




FIG. 1: A section of the pyrochlore lattice, showing its 
two-sublattice tetrahedral structure. At zero-magnetization, 
tetragonal lattice distortion favors the spin configurations of 
the type shown in (a), while at half-magnetization trigonal 
lattice distortion favors the configuration (b). AF bonds are 
marked with black and FM bonds with red lines. 
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HgCr 2 4 EH. 

We take as a starting point the Hamiltonian 



H = 



(id) 
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where the summation runs over the nearest neigh- 
bor bonds of a pyrochlore lattice and pij is the change 
in distance between neighboring spins S j and S j , relative 
to the equilibrium lattice constant. We assume the ex- 
istence of a linear regime in which exchange interactions 
and elastic energies depend only on the distance between 
lattice sites. It is convenient to introduce a single dimen- 
sionless parameter 6 = Ja 2 /K to measure the strength 
of the spin-lattice coupling. As written, the antiferro- 
magnetic exchange interaction J, elastic constant K and 
a spin-lattice coupling a, are all taken to be positive. 

The pyrochlore is a bipartite network of corner sharing 
tetrahedra. This means that, if we neglect coupling to the 
lattice, we can write the energy per spin as 



W = 4J]T (M-— ^ 



tctr. 



8 J / 



h 2 
16J 



const. 



(2) 



where the sum runs over all the tetrahedra and 
M = (Si + S 2 + S3 + S4)/4 is the magnetization per 
site. Clearly, the energy can be minimized by ensuring 
M = h/8 in each tetrahedron. However this local con- 
straint does not select any one ground state, but rather 
a vast manifold of states. If we consider classical spins, 
one angular variable remains undetermined per tetrahe- 
dron, and the magnetization is linear in h right up the 
saturation field of h = 8 J. 

Coupling to the lattice provides a very efficient mech- 
anism for lifting this degeneracy. Since bond energies 
vary linearly with p^j, while elastic energies increase as 
pf j , at any given value of magnetic field the system can 
always gain energy by ordering the spins and distorting 
the lattice. In this sense the Hamiltonian can be 
thought of as a three dimensional generalization of the 
spin-Peierls problem. The system gains the most en- 
ergy from distorting bonds for which SjSj takes on its 
extremal values. For this reason, coupling to the lat- 
tice tends to favor collinear spin configurations and, for 
quantum spins, bond singlets (see e.g. 0]). Our goal is 
to understand which states emerge from this competition 
between applied field and frustrated AF interactions, and 
for what range of fields they are stable. 

For the sake of simplicity, and in the spirit of earlier 
work 0, 0], we shall restrict our analysis to uniform 
spin and lattice order with crystal momentum q = 0. It 
is instructive to further simplify the problem by treating 
the spins as classical vectors, in which case we can safely 
neglect all states which are odd under the inversion It 
which exchanges the two tetrahedron sublattices. Both 
of these approximations can be relaxed at will. 




FIG. 2: Phase diagram as a function of magnetic field h and 
dimensionless coupling constant b. Solid lines denote first and 
dashed lines second order transitions. Spin configurations and 
irrep of the order parameter in each phase is also shown. 



Under these assumptions, the system must have four 
sublattice order, and we can find the ground state of 
Eq. Q by a straightforward minimization of energy with 
respect to bond length. This is equivalent to sol ving a 
Heisenberg model with biquadratic-bilinear terms [l5j : 



(3) 



(id) 



Our results are summarized in the phase diagram Fig. [21 
with the corresponding magnetization curves shown in 
Fig. 01 For small h, the lattice has overall tetragonal sym- 
metry, with tetrahedra distorted so as to have two long 
ferromagnetic (FM) and four short canted AF bonds. 
This is broadly compatible with the experimentally ob- 
served ground state of ZnCr 2 04 In this regime the 
magnetization of the system remains linear, but with re- 
duced slope M w ft/(8J(l + 26)). 

For h w 3J, and b > 0.05 the system makes a first 
order transition into state with exactly 3 up and 1 down 
spins per tetrahedron, i.e. M = 5/2, regardless of h. 
In this half-magnetization "plateau" phase, each tetra- 
hedron has three long FM and three short AF bonds, 
giving rise to an overall trigonal lattice distortion. For 
any finite value of b the plateau is extremely broad. Its 
width shrinks linearly as b — > 0; for hj J < 4 — 166 we find 
a transition into a coplanar 2:1:1 canted state with mixed 
E and T 2 symmetry, and for hj J = 4 + 86 a transition 
into a coplanar 3:1 canted phase with trigonal symmetry 
for h/J = 4 + 8b. Both transitions are of second order. 

Finally, for 8J > h > 6J, there is a transition into 
a fully saturated FM in which the lattice regains overall 
cubic symmetry. In the absence of longer range exchange 
interactions, the two lowest lying spin wave branches of 
the FM phase are local in character and dispersionless. 
For 6 < 3/38, the transition from the 3:1 canted phase 
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FIG. 3: Magnetization as a function of magnetic field for 6 = 
(straight line) to b — 0.24 (robust plateau) in steps of 0.03. 



into the FM is of second order, and occurs on the line 
h/J = 8 — 166. For 6 > 3/38 the transition becomes 
first order, and for 6 > 1/4 it is energetically favorable 
to make a transition directly from the collinear plateau 
phase into the FM. 

In order to understand why these particular phases are 
stable, we now turn to the symmetry analysis. For clas- 
sical spins, coupled to uniform lattice distortion, we need 
only consider the symmetries 7^ of a single tetrahedron 
|l6j. The bond variables pij , which describe changes in 
the length of the six edges of the tetrahedron, transform 
according to the Ai, E, and T2 irreducible representations 
(irreps) of T d : 
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{pe,i,Pe,2} irrep (which includes tetragonal 
distortions of the lattice), opposing pairs of bonds deform 
with the same sense. In the p Tl — {pr 2 .i, Pt 2 ,2, PT2.3} ir- 
rep (which includes trigonal distortions ), opposing pairs 
of bonds deform with the opposite sense. Exactly anal- 
ogous representations for spins can be obtained by sub- 
stituting with S.;Sj for pi j and A Al for p Al , etc. 

In terms of these variables, the Hamiltonian for a single 
tetrahedron embedded in the lattice is given by 



In the p E 



H = 2V6JA A 
+K (pi H 



2a J (A A p A + A E p E 
p\ 2 ) - 4hM , 



At 2 Pt, 



(•5) 



where Arp r and = p R p R are second order invariants 
of the R = Ai, E, and T2 irreps. An analysis of Eq. © in 
the absence of magnetic field was given in Refs. [l3| and 
[14), In all of the cases considered by these authors, only 




FIG. 4: Maximal values of the second order invariants Aj| 
and Ay 2 for classical spins as a function of the magnetization 
per site M. A ridge (dashed line) divides three dimensional 
from coplanar spin configurations. Spin ordering patterns are 
shown for the symmetrical cases of Al = and K\ 2 = 0. 



the E irrep is relevant. However, once a magnetic field 
is applied, both the Ai and T2 irreps have an important 
role to play. For classical spins, Eq. JSJ reduces to 

E a = 2 j(V6A Al - 6 Al A Ai - b E A 2 E - 6y 2 A T2 ) - 4hM . (6) 

For pure nearest neighbour interaction [c.f. Eq. JQl] 
^Ai = = bj 2 = b. In general, however, these parame- 
ters need not be equal. 

Using the fact that A Al = 8 (AI 2 - 1 /4)/\/6, we see that 
the lowest energy configuration at a given magnetization 
is that for which 6eA e + 6j 2 A T2 takes on its maximal 
value. The surface of maximal values of these second 
order invariants is shown in Fig. 01 The limiting cases 
Aj 2 — > or Ae — > have a simple analytic form 



max Ap 

A? =0 E2 



max At 



16(1-M 2 )73 



(1 - M 2 



32 



if 1/2 < M < 1 
if < M < 1/2 



(7) 



(8) 



and the stability of the half-magnetization plateau orig- 
inates in the sharp cusp in the maximal value of A 2 - as 
a function of M. Provided that 6e < 26y 2 , this singular- 
point (which corresponds to the trigonal lattice distor- 
tion shown in Fig.^ 3 ), minimizes the energy for a finite 
range of values of magnetic field, and the phase diagram 
is qualitatively that of Fig. [21 Conversely, for values of 
6e > 26y 2 , T2 order is not realized for any value of h. 
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As the magnetization of the system increases, so will 
the average bond length, and for classical spins coupled 
to a uniform lattice distortion, the volume of the unit cell 
is a monotonically increasing function of applied field 

SV [3 1 2 7 , 3 2 

1/ = V2 PAl + 2^"4 PE "4 PT2 - (9) 

From this expression we learn that: (i) as pAi oc Aa i; 
the jumps in magnetization will have their counterpart 
in volume change; (ii) the application of hydrostatic pres- 
sure will soften the E mode relative to the T2 mode, thus 
extending the region of the tetragonal phase. 

The scenario which we have presented is the simplest 
under which lattice distortion can stabilize a magnetiza- 
tion plateau in a spinel oxide. Needless to say, the situ- 
ation in a real spinel oxide may be much more complex 
than that described above. Naive estimates suggest that 
there are of order 1.3 N / 2 collinear states with 3 up and 
1 down spins in each tetrahedron, where N is the num- 
ber of spins in the lattice. In principle, any of these may 
couple to phonon modes. Arbitrarily complex exchange 
interactions, competing with arbitrarily complex elastic 
energies, may give rise to an arbitrarily complex plateau 
state — or none at all. None the less, our model captures 
the essential features of the high field magnetization of 
CdCr 2 04 and HgCr 2 04 0. For these Cr spinels t 2g 
shell is full and the e g shell empty, so we may safely ne- 
glect the effects of orbital degeneracy 0. The theory 
may also be applicable to ZnCr204, although the high 
value of J in this compound makes verification difficult. 

Chromium spinels exist as chalcogenides as well as 
oxides 0. The chalcogenides have weaker and predom- 
inately FM interactions, and are therefore less likely a 
priori to exhibit a half-magnetization plateau. However 
the competition between FM and AF interactions is al- 
ready present in oxide materials. We have examined the 
role of next-nearest neighbour exchange J2 within spin 
wave theory, the minimization of the energy of small clus- 
ters, and an extension of the symmetry analysis presented 
above. For AF J2, bj 2 > &e, and the stability of the trig- 
onal half-magnetization plateau phase is enhanced. FM 
J2, on the other hand, drives the system towards lattice 
distortions (and associated plateau states) at finite q. It 
is also worth noting that the strength and sign of ex- 
change interactions in oxides and chalcogenides can be 
very sensitive to bond angle. As written, the Eq. (JTJ 
does not allow for the rotation of neighboring tetrahedra 
at fixed bond length. Such modes will be important at fi- 
nite q, and may lead to a magnetostriction (i.e. decrease 
in the volume of the unit cell at in applied field). The 
necessary extension of our theory to treat these cases is 
in principle possible, but lies beyond the scope of the 
present letter. 

We conclude with a few comments on the role of fluctu- 
ations. In constructing a theory of a half-magnetization 



plateau stabilized by lattice distortion we have assumed 
static spin order. Since both quantum and thermal fluc- 
tuations in frustrated magnets favour collinear spin con- 
figurations j these will further contribute to the sta- 
bility of the magnetization plateau. If consider Eq. as 
an effective hamiltonian, the coefficient b will have con- 
tributions fluctuation effects as well as lattice distortion 
[2JJ] . It may also have contributions of a purely electronic 
origin 21]. However, we have performed monte carlo 
simulations of Eq. (JJJ for classical spins in the absence 
of coupling to the lattice, and these suggest that any 
plateau stabilized by "order from disorder" effects alone 
will be at least as fragile as that seen in the Kagome lat- 
tice || . The remarkable width of the half-magnetization 
plateaus observed in CdCr 2 04 and HgCr 2 04 — which ex- 
tend over many tesla — leads us to the conclusion that 
lattice distortion plays a crucial role in these systems. 
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